]> permondes.de Git - Analog_Engine.git/blob - scripts/TP1 07.19 Falling Mass.SALSA
Applications: Longitudinal oscillating chains; PDEs
[Analog_Engine.git] / scripts / TP1 07.19 Falling Mass.SALSA
1 # Wiring mache ich von links nach rechts !!!
2
3 ...IDENTIFICATION DIVISION
4 ...PROGRAM-ID FallingParticle
5 ...VERSION 20240225
6 ...COMMENT Calculate the deflection from the vertical caused by the Earth's rotation of a particle falling freely from rest from a height h.
7 ...COMMENT Differential Equations:
8 ...COMMENT x''=-bz'+ay' # x-axis is along latitude, directed to east
9 ...COMMENT y''=-ax' # y-axis is along longitude, directed to north
10 ...COMMENT z''=-g+ax' # z-axis is perpendicular to the surface of earth
11 ...COMMENT g: gravitational acceleration = 9,81 m/s²
12 ...COMMENT a: 2*omega*sin(phi)
13 ...COMMENT b: 2*omega*cos(phi)
14 ...COMMENT omega: rotation velocity of the earth = 2*pi/day
15 ...COMMENT phi: Latitude of location (0-90°)
16 ...COMMENT Initial Condition: z(0)=h
17 ...COMMENT The full solution requires 6 INTEGRATORs, Anabrid-THAT just has 5. The deflection to longitude (y) is neglegible and can be omitted (marked #*).
18 ...COMMENT It could also be solved in a separated algorithm omitting x.
19
20 ...ENVIRONMENT DIVISION
21 ...TIMEBASE 1ms
22
23 ...DATA DIVISION
24 SET COEFFICIENT.1 TO AY # 2*2pi/day*sin(phi)
25 SET COEFFICIENT.2 TO AX # = AY
26 SET COEFFICIENT.3 TO B # 2*2pi/day*cos(phi)
27 SET COEFFICIENT.4 TO G # gravitational acceleration = 9,81 m/s²
28 SET COEFFICIENT.5 TO H # height h
29 SET OUTPUT.X TO x
30 SET OUTPUT.Y TO y
31 SET OUTPUT.Z TO z
32
33 INITIALIZE H by -1 TO -h # same as COMPUTE -1 TIMES H TO -h
34 INITIALIZE G by +1 TO g
35
36 INTEGRATE 1*-bz',1*ay' TO -x' # Input is x''
37 INTEGRATE -x' TO x
38
39 INTEGRATE y'' TO -y'
40 INTEGRATE -y' TO y
41
42 INTEGRATE z'' TO -z'
43 INTEGRATE -z', IC:-h, LIMIT:(z >= 0) TO z
44
45 -x' * AX -> -ax'
46 -ax' = y''
47 -y' * AY -> -ay'
48 INVERT -ay' TO ay'
49 -z' * B -> -bz'
50 ADD 1*-ax', 1*g TO -g+ax'
51 -g+ax' = z''
52
53 ...OPERATION DIVISION
54 ...MODE REPEAT
55 ...OP-TIME 7,3ms